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ABSTRACT 

Three planets have been directly imaged around the young star HR 8799. The planets are 5-13 Mjup and 
orbit the star at projected separations of 24-68 AU. While the initial detection occurred in 2007, two of the 
planets were recovered in a re-analysis of data obtained in 2004. Here we present a detection of the furthest 
planet of that system, HR 8799 b, in archival HST/NICMOS data from 1998. The detection was made using 
the locally-optimized combination of images algorithm to construct, from a large set of HST/NICMOS images 
of different stars taken from the archive, an optimized reference point-spread function image used to subtract 
the light of the primary star from the images of HR 8799. This new approach improves the sensitivity to planets 
at small separations by a factor of ~10 compared to traditional roll deconvolution. The new detection provides 
an astrometry point 10 years before the most recent observations, and is consistent with a Keplerian circular 
orbit with a ~70 AU and low orbital inclination. The new photometry point, in the F160W filter, is in good 
agreement with an atmosphere model with intermediate clouds and vertical stratification, and thus suggests the 
presence of significant water absorption in the planet's atmosphere. The success of the new approach used here 
highlights a path for the search and characterization of exoplanets with future space telescopes, such as the 
James Webb Space Telescope or a Terrestrial Planet Finder. 
Subject headings: planetary systems — techniques: image processing 



1. INTRODUCTION 

After more than a decade marked by the great success of 
the radial velocity and transit planet detection techniques, 
the direct imaging technique has finally made its grand en- 
try on the scene in la te 2008 thro ugh a se ries of exciting 
exoplanets discoveries dMarois et a l. 2008; Kalas et al.ll2008l; 
iLafreniere et al.l l2008t iLagrange e t al.i i2008i) . including the 
spectacular discovery of the multiple-planet system HR 8799. 
This system showcases three planets of mass 5-13 Mjup or- 
biting at projected separations of ^^24, 38 and 68 AU. It is 
located 39.4 pc away from the Sun and appears to be seen 
nearly face-on. The planets all orbit the star in the same di- 
rection and likely in the same plane, consistent with forma- 
tion within a circumstellar disk. The host is a 30-160 Myr- 
old star of spectral type A5, mass ^1.5 Mq and luminosity 
~4.9 Lq; it is also classified as both 7 Dor and A Boo types. 
A large infrared excess is detected at ~100 /im, suggesting 
that the planetary system is surrounded by a large dust disk. 
The reader is referred to IMarois et al.l (l2008b and references 
therein for further details on this system. 

As the first multiple-planet system directly imaged, 
HR 8799 offers many new possibilities to advance our un- 
derstanding of planets. First, multi-wavelength photometry 
and spectroscopy of these three planets, which are almost cer- 
tainly of the same age and metallicity, will be of great value 
for the validation and calibration of both evolutionary and at- 
mospheric models of giant planets. Also, with good measure- 
ments of the orbits of the planets, it could be possible to inde- 
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pendently constrain the ir masses through dynamical stu dies, 
as already attempted by lFabrycky & Murray-Clayl(l2008i) . 

Besides its high scientific interest, the discovery of the 
HR 8799 system represents a great technical achievement 
given the brightness ratios and angular separations between 
the planets and the star. The main obstacle to direct imaging 
of exoplanets, for both space- and ground-based telescopes, 
is the scattering of stellar light by irregularities on optical 
surfaces which creates a halo of bright quasi-static speckles 
around the stellar core, masking the underlying fainter plan- 
ets. The planets around HR 8799 were initially discovered us- 
ing adaptiv e optics and the ang ular differential imaging (ADI) 
technique (iMarois et al.ll2006l) . which uses the natural field- 
of-view rotation of an alt-az telescope to discriminate plan- 
ets from quasi-static speckles. Effectively, this techniques al- 
lows to construct a high-fidelity image of the stellar point- 
spread function (PSF) that does not contain the signal of any 
eventual planets; this reference PSF image is subtracted from 
the target image to remove the light from the star while pre- 
serving that of any eventual planet. This technique, coupled 
with the locally-optimized combination of images (LOCI) al- 
gorithm (ILafreniere et al.ll2007b ). has proved to be the most 
efficient approach to de tect planets from the ground (see e.g. 
ILafreniere et"ai]|2007ah . 

Following the initial detection of HR 8799 b and c in Oc- 
tober 2007, a careful re-analysis of adaptive optics observa- 
tions of HR 8799 obtained in July 2004, using improved PSF 
subtraction algorithms, allowed to recover both planets, thus 
providing a baseline of four years relative to the most re- 
cent observations of September 2008. This four-year baseline 
firmly establishes common proper-motion of the two outer- 
most planets and clearly reveals their orbital motion. How- 
ever, this baseline is still quite small compared to their orbital 
periods (P > 200 yr), and astrometric measurements over a 
longer baseline would be highly valuable in better constrain- 
ing their orbits, and consequently the star and planets masses. 
This will require several years of observations, but could be 
achieved more rapidly should a similar re-analysis be applied 
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successfully to even earlier data. 

HR 8799 was observed with HST/NICMOS in 1998 as part 
of a direct imaging survey for massive planets around young 
nearby stars (program 7226, PI Eric Becklin); the main re - 
sults of this survey were published in lLowrance et alJ (l2005h . 
Although these observations have, in principle, sufficient an- 
gular resolution and flux sensitivity to see the planets of the 
HR 8799 system, the scattered light from the bright primary 
prevented their actual detection. This scattered light can be 
partly removed by subtracting images obtained at two dif- 
ferent spac ecraft roll angles, the so-called roll-deconvolution 
technique (ISchneider & Silverstonell2b03h . but this is insuffi- 
cient to reveal the planets because of important PSF evolution 
between the two roll angles. As part of the same HST pro- 
gram, many other stars were observed using the same instru- 
mental configuration as for HR 8799, and as a result several 
images d isplay PSFs very similar to that of HR 8799. As sug- 
gested in lLafreniere et al.l (l2007bh . such observations provide 
a very interesting basis for constructing an optimized refer- 
ence PSF image us ing the LOCI a lgorithm. The LOCI al- 
gorithm, detailed in lLafreniere et a l. (2007b), determines the 
coefficients needed to linearly combine several reference im- 
ages into an optimized reference PSF image whose subtrac- 
tion from the target image will minimize the residual noise. 
An important and powerful feature of this algorithm is its 
flexibility to optimize the PSF subtraction locally over sev- 
eral sub-sections of the image. Applied to the HST data set 
mentioned above, the LOCI algorithm could potentially re- 
duce the scattered light of HR 8799 sufficiently to reveal the 
planets. With this in mind we have re-analyzed these data 
and here we report the successful detection of the furthest of 
the three planets, HR 8799 b, thus extending astrometric mea- 
surements of this planet to ten years. Incidentally, this means 
that an exoplanet could have been discovered by direct imag- 
ing only a few years after the discovery of 51 Peg b by the 
radial velocity technique. 

2. DATA SET AND ANALYSIS 

All the data used in this study come from HST program 
7226, which was a survey for giant planets and brown dwarfs 
around young nearby stars. The observations used the NIC- 
MOS medium resolution camera NIC2, the F160W filter and 
the 0.63" diameter occulting spot to reduce the diffracted light 
from the primary star. Typically three images of each target 
were obtained at each of two roll angles differing by 29.9°. 
A total of 39 targets were observed for this program. The 
pipeline-reduced images used here were retrieved from the 
HST archive housed at the CADC.^ 

As mentioned above, the LOCI algorithm was used with 
this large set of well-correlated images to construct and sub- 
tract an optimized PSF for each of the six images of HR 8799. 
A priori, the set of reference images used by the LOCI algo- 
rithm could include all images of targets that are not HR 8799, 
but in practice we have omitted some targets from the set be- 
cause they showed a second point source within a separation 
of 3" from the star, were significantly fainter than HR 8799, or 
were not centered behind the coronagraph. The final set of ref- 
erence images included 203 images from 23 unique sources. 
When analyzing a given image of HR 8799, we have added 
to this set the three images of HR 8799 obtained at the other 
roll angle. Before execution of the LOCI algorithm, all ref- 
erence images were registered to the target image based on a 

^ |http : //www2 ■ cadc-ccda ■ hia-iha ■ nrc-cnrc ■ gc ■ ca/cadc/ 




Fig. 1. — (a) Image of HR 8799 at first roll angle, (b) Image at the first 
roll angle minus image at the second roll angle, (c and d) Residual images, 
at both roll angles, after subtraction of an optimized reference PSF image 
obtained using the LOCI algorithm, see text for more details. The display is 
from to 100 counts s"' for panel a and from -0.5 to 2.5 counts s"' for panels 
h-d. Each panel is 4.8" on a side. The dotted circles in the bottom two panels 
mark the inner and outer limits of the optimization annulus used for the LOCI 
algorithm. The two short radial lines enclose the position of HR 8799 b. 
To produce the images shown in panels c and d, a second optimization was 
performed for the half annulus not containing HR 8799 b, such that the full 
annulus is representative of the residual noise. 

cross-correlation of the secondary mirror support diffraction 
spikes. The algorithm was first applied with the specific goal 
of detecting HR 8799 b, and thus a single optimization region 
was used. This optimization region was defined as a half an- 
nulus extending radially from 1.275" to 2.175" and centered 
on, but excluding pixels within a radius of 5 pixels from the 
calculated position of HR 8799 b. The exclusion was used 
to prevent the algorithm from trying to subtract the compan- 
ion itself, thus biasing the determination of the coefficients 
of the linear combination and affecting the residual flux of 
the companion. The area of the optimization region, ~900 
pixels, is equivalent to 225 PSF cores (i.e. the characteristic 
size of a speckle). This number of "degrees of freedom" in 
the optimization region may be compared with the number of 
different variables, 46 different observations (23 sources each 
at two roll angles), in the linear combination of the optimum 
reference image. The different observations do not represent 
free parameters of the model PSF, however, as they are highly 
correlated with each other 

The companion HR 8799 b was detected in all six resid- 
ual images and its position between the two rolls angles is 
precisely equal to the spacecraft rotation angle applied. The 
three residual images at each roll angle were co-added and the 
results are shown in Fig.[T] At the peak pixel of the compan- 
ion in these co-added images, the detections are at the levels 
of lOcr and 6<t for the first and second roll angles, respec- 
tively. We have verified that the residual noise in the opti- 
mization region closely follows a Gaussian distribution, thus 
the significance of the detection is high. The improvement in 
PSF subtraction provided by the LOCI algorithm is obvious 
from Fig. [T] in which the much higher residuals of the roll- 
subtracted image is apparent. At the separation of HR 8799 b, 
the residual noise in the LOCI-subtracted image is nine times 
lower than for the classical roll-deconvolution image (panel b 
of Fig. [U. It is interesting to note that the LOCI subtraction 
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TABLE 1 
HR 8799 b astrometric measurements 



Epoch 


Aa (mas) 


AS (mas) 


Sep (mas) 


PA(°) 


1998.8285 


1411±9 


986 ±9 


1721 ±12 


55.1 ±0.4 


2004.5337 


1471 ±5 


884 ±5 


1716±7 


59.00 ±0.23 


2008.6267 


1527 ±2 


800 ±2 


1724 ±3 


62.35 ±0.09 



Note. — Ep. 1998 from this work; others from lMarois et alj 120081) . 



almost reaches the PSF photon noise. At the separation of 
HR 8799 b, we estimate that the total residual noise is ^^1.8 
times more important than the local photon noise. 

To confirm that the above detection is not an artifact from 
the image processing done, we have repeated the analysis 1) 
without excluding the region containing the companion from 
the optimization half-annulus, and 2) using only half of the 
available reference images. The companion was still well de- 
tected in both cases, although at lower signal-to-noise ratios 
(S/Ns) as expected. We have repeated the above analysis with 
an optimization region designed specifically for HR 8799 c, 
at ^0.9", but the results were inconclusive. The same applies 
for d, which is at an even smaller separation. 

3. RESULTS 

The astrometric and photometric analysis of the companion 
was done using mo del PSFs gen erated with the TinyTim^ soft- 
ware (version 7.0) (lKrisllll993h . Model PSFs appropriate for 
the NIC2 pre-cryocooler camera with the F160W filter were 
generated for each of three approximate positions on the de- 
tector: the central star, the companion at the first roll angle, 
and the companion at the second roll angle. All model PSFs 
were generated for a width of 10" and a pixel oversampling 
factor of 9 to minimize interpolation effects when shifting im- 
ages. Accordingly, during the analysis all shifts were done on 
the oversampled images, which were then binned 9x9 pixels 
to match the actual NIC2 plate scale. Initially, the model PSFs 
are normalized to a total flux of one. 

For each image, the appropriate companion model PSF was 
first spatially shifted and intensity scaled over a grid in dx, 
dy, and flux. Then each model PSF in this 3D parameter 
space was subtracted from the image and the residual noise in 
a 7 X 7-pixel box centered on the companion was computed. 
The combination of dx, dy, and flux that yielded the minimum 
residual noise was found and provided the coordinates and 
flux of the companion. The uncertainty on the position of the 
companion, 5 mas, was estimated from the dispersion of the 
measurements made in all individual images. The center of 
the star was determined by shifting the appropriate model PSF 
to maximize the cross-correlation of its diffraction spikes with 
those of the actual image, as well as by a cross-correlation of 
the observed PSF diffraction spikes with themselves after a 
rotation of 180° about the inferred center; the maximum dif- 
ference between these two centroid measurements, ^0. 1 pixel 
or 7.5 mas, was taken as an estimate of the centroid accuracy. 
The pixel coordinates of the star and companion were con- 
verted to RA and DEC using the astrometry solution defined 
in the image FITS headers. A S/N-weighed mean of the rel- 
ative position of HR 8799 b over the six images was finally 
calculated; the results are indicated in Table[T] 

Since the model PSFs span 10", the companion flux found 
by the above procedure is effectively the same as that in 
an infinite aperture. The information in the NICMOS Data 

^ |http : / /vivivi ■ St sci ■ edu/sof tware/tinytim/tinytim. html| 
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Fig. 2. — Measured relative positions (with associated en'ors) of HR 8799 b 
at the three available epochs. The solid curve shows a Keplerian circular 
orbit with a = 71.4 AU, i = 19.2°, and position angle of ascending node Q = 
131.3°; the open diamonds mark the corresponding planet position at the 
three epochs. The for this orbit is 3.19. 



TABLE 2 
HR 8799 b photometric measurements 



Filter 




Magnitude 




Measured 


Intermediate clouds 


Extreme clouds 


/ 


19.28±0.16 


19.42 


19.66 


F160W 


18.54±0.12 


18.58 


18.31 


H 


17.85 ±0.17 


18.30 


18.08 




17.03 ±0.08 


17.43 


16.91 


L' 


15.64 ±0.11 


15.24 


15.59 



Note. — F160W from this work; others from lMarois et alj <2008l) . 



Handbook v7.0 was then used to convert this flux (20.5 ± 
2.5 counts s"') to physical units, yielding (4.2±0.5) x 10"^ Jy, 
or a magnitude of 18.54 ±0.12. This magnitude is also in- 
dicated in Table [2] along with the measurements reported by 
iMarois et al.l (12008,) and the expected magnitudes for two dif- 
ferent atmosphere models (see ^for more detail on models). 

4. DISCUSSION 

As visible in Fig. |2] the new position is consistent with the 
previous ones, and still suggests a nearly circular orbit seen 
close to face-on. We have done very simple orbital fits to 
verify that the data are consistent with true Keplerian orbits. 
For simplicity, we have considered only circular orbits and as- 
sumed that the stellar mass is precisely equal to 1.5 Mq, and 
then we explored a range of semimajor axes (a = 60-100 AU) 
and inclinations (/ = 0^5°). The best fits, with s. ^ 3.2, 
were found for a ^ 68-74 AU and / ^ 13-23°; an example 
of such fit is shown in Fig.|2] This range of inclination is in 
line with what is expected for the equato rial plane of the star 
based on its measured v sin/ (37.5 km s-i . iGrav & KavJT999l) 
which, for the range of true rotation velocity of A5 stars (100- 
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Fig. 3. — Synthetic spectrum from a 820 K, logg = 4.0 model atmosphere 
with veilically stratified clouds. Filled symbols are the band-integrated fluxes 
and observations are shown as open symbols with 3cr error bars. Filter re- 
sponse curves are plotted below the spectrum. 

300 km s-i, iRover et al.ll2007h . would yield / = T-IT. We 
refrain from carrying out more detailed orbital fits as the lim- 
ited time baseline and astrometric precision available would 
prevent us from reaching useful constraints on the orbits. 

The low luminosity and red near-infrared color of 
HR 8799 h{J -K^ = 2.25) are indicative of atmospheric dust, 
but constraining the dust properties (composition and distri- 
bution) requires a detailed comparison of its SED with var- 
ious model atmosphere predictions. Figure [3] compares the 
new F160W photometric point, along with earlier 7, H, 
and L'-band photometry, to an intermediate (vertically strat- 
ified) cloud model (Barman et al., in prep.) with parame- 
ters selected by comparing the planet's age and luminosity to 
substellar evolution tracks (iMarois et al.ll2008l) . The observed 
and model photometry are listed in Table |2] The F160W fil- 
ter is wide enough to incorporate a substantial fraction of the 
water band between the J and H bands. Since the depths of 
the water bands are greatly reduced as dust content increases, 
this filter, and other NICMOS band-passes, have greater po- 
tential for constraining the atmospheric dust content than the 
standard ground-based near-IR windows. As can be seen in 
Fig. [3] (and Table |2]i, the new F160W photometric point is 
in excellent agreement with this intermediate cloud model, 
and about 2a fainter than an extreme dusty atmosphere hav- 
ing clouds that blanket most of the atmosphere. While the 
overall agreement with the extreme cloud atmosphere model 
appears to be relatively good, based on Table |2] this model 
requires Tgff = 1600 K (cf '-^SOO K for the intermediate cloud 
model) and a very small radius (^4 Earth radii) to match the 
observed luminosity. Both of which are inconsistent with for- 
mation and evolution models. 



5. CONCLUDING REMARKS 

The work presented here demonstrates that the LOCI algo- 
rithm can take HST data a step further, extending the sensi- 
tivity to planets beyond what has been achieved before. The 
good performance of the LOCI algorithm with HST data is 
due in large part to the high stability of the optical aberra- 
tions of HST and the consistent, accurate pointing between 
the different visits. Over the years, a large archive of HST ob- 
servations aimed at detecting exoplanets has been assembled. 
Based on the present results, it would be interesting to look 
again at this set of data, using the LOCI algorithm, to see if 
any other sources have been missed by previous searches. On 
a related note, it could be interesting to attempt a re-analysis 
of archival ground-based AO observations using an approach 
similar to that used here. However, this could turn out be 
much less efficient than for HST given the larger evolution of 
the PSF structure for ground-based telescopes. 

Looking into the future, the approach presented here should 
definitely be an integral part of the strategy for the search 
and characterization of planets with future space telescopes, 
such as the James Webb Space Telescope (JWST) or a Ter- 
restrial Planet Finder JWST, for instance, is expected to be 
relatively stable in temperature and should not suffer from 
the breathing problem experienced by HST. Its PSF should 
therefore be more stable than that of HST and the LOCI al- 
gorithm will likely perform extremely well, especially for ob- 
servations obtained within a giv en wavefront adjust ment cam- 
paign, typically every 14 days ([Gardner et alJ 12006 ). A pre- 
liminary an alysis of the JWST PSF temporal evolution done 
by Makido n et al.l(l2008l) indeed suggests that PSF subtraction 
should perform well. For this approach to work, enough atten- 
tion should be paid to the accurate positioning of the different 
targets, in particular for observations made with the corona- 
graphic occulting masks. Having an efficient PSF subtraction 
strategy will be absolutely necessary for JWST given its com- 
plicated, speckled PSF structure arising from its segmented 
pupil, in addition to the speckles arising from the unavoid- 
able optical aberrations. While roll-deconvolution could per- 
form well for JWST in the absence of telescope breathing, the 
LOCI approach remains very interesting given the observa- 
tory's maximum roll angle of only ±5°, which will severely 
limit the use of roll-deconvolution at small separations. 
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